Anomalous thickness dependence of the Hall effect in ultrathin Pb layers on Si(lll) 
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The magnetoconductive properties of ultrathin Pb films deposited on Si(f f f) are measured and 
compared with density-functional electronic band-structure calculations on two-dimensional, free- 
standing, 1 to 8 monolayers thick Pb(lll) slabs. A description with free-standing slabs is possible 
because it turned out that the Ifall coefficient is independent of the substrate and of the crystalline 
order in the film. We show that the oscillations in sign of the Hall coefficient observed as a function 
of film thickness can be explained directly from the thickness dependent variations of the electronic 
bandstructure at the Fermi energy. 
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The electronic properties of crystalline bulk materials 
are well understood. However, as soon as the lateral 
dimension in any direction is reduced to a few atomic 
distances, each electron energy band splits into a set of 
discrete subbands and the spatial confinement leads to 
quantum-size effects (QSE). The quantum-size confine- 
ment affects in particular the electron transport prop- 
erties. In an early theoretical approach Trivedi and 
AshcroftE considered transport of electrons in a metal 
film, confined in a two-dimensional potential well and 
scattered by impurities in the film and/or by surface 
roughness. They predicted a saw-tooth like variation of 
the conductivity with the film thickness with a period 
of half the Fermi wavelength. Such oscillations have not 
been observed in the experiment so far. Another source 
of QSE-induced conductivity oscillations is the layer-by- 
layer growth of, e.g., Pb films on Si(lll)(7x7) or Si(lll)- 
Au(6 X 6) when the in-plane conductivity shows oscilla- 
tions with the period of the interlayer spacing, associ- 
ated with the degree of roughness of the surfaceE'0 An 
even more spectacular, but so far unexplained, observa- 
tion was the reversal of sign in the Hall coefficient Rh 
of epitaxial Pb films on Si(lll), grown on a surface that 
was pre-covered with about one monolayer (ML) of Ag 
which shows a (Vs x \/2>) reconstruction.a 

In this report, we address the general issue of charge 
transport and Hall effect in ultrathin metallic layers both 
theoretically and experimentally. We focus on Pb films 
deposited on Si(lll), although the conclusions should 
be generally valid for ultrathin metallic films in the re- 
gion of QSE. In order to clarify the underlying mech- 
anisms, we combined in our study theoretical investi- 
gations by density-functional calculations with measure- 
ments of conductivity and of magnetoconduction of ultra- 
thin Pb films on a (7x 7) reconstructed Si(lll) substrate. 

Pb films with a thickness between 1 and 12 monolayers 
(ML) have been deposited in UHV on a well annealed 
Si(lll)— (7 X 7) substrate at low temperatures (15 to 
25 K). The thickness, crystalline order and defects have 
been monitored with higbpresolution low-energy electron 
diffraction (SPA-LEED).B 4 Mo contacts have been pre- 



dcposited for the measurements. After the pseudomor- 
phic first Pb layer, up to 4 Pb layers were amorphous, 
whereas thicker layers always turned out to be crystalline. 
By varying temperature and thus the growth conditions, 
the influence of order on the conductive and magnetocon- 
ductive properties was tested. For details on film prepa- 
ration and ordering in the films, see Ref. ^ 

Before starting the Hall effect measurement the films 
have been annealed close to temperature where the de- 
crease of conductance indicated a break-up of the film 
into non-contiguous islands. Even after annealing, films 
with a thickness d < 4ML were strongly disordered 
whereas thicker films grew epitaxially with the lattice 
constant of bulk Pb. While the conductance measure- 
ments, which provide information on the charge-carrier 
scattering, were to some extent sensitive to the degree 
of order, the oscillatory behavior of the Hall coefficient 
turned out to be remarkably insensitive. Therefore, it 
seems to be justified as a first approach to compare the 
experimental results of the Hall coefficient with those ob- 
tained theoretically for free-standing Pb slabs. Indeed, as 
we will show below, such calculations are able to identify 
the basic origin of the oscillations. 

Conductance and magnetoconductance measurements 
were performed in a DC magnetic field up to 4 T at 
various temperatures down to about, 7 K. In the van-der- 
Pauw arrangement for conductance^ it has been checked 
that the magnetoconductance was completely symmet- 
ric. In the contact arrangement for Hall effect there was 
a sizable contribution from magnetoconductance due to 
deviations from a perfect quadratic arrangements of the 
contacts. Therefore the asymmetry was taken as the con- 
tribution of the Hall effect to the measured voltage. The 
Hall coefficient has been derived from this asymmetry at 
±4T. 

In the model calculations we treat the Pb layers as free- 
standing slabs and neglect, as mentioned, the effect of the 
substrate. Within the model we calculate, as a function 
of layer thickness, the densities of states at the Fermi en- 
ergies, Fermi velocities, and effective-mass tensors. From 
these quantities, together with the experimental data for 



the conductivity, cr™^, the elastic scattering times are 
calculated and the Hall coefficient, Rh is derived. 
The Hall coefficient in the low-field limit is 
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where ctq and an a-i'e the electrical and Hall conduc- 
tivities of a two-dimensional slab. Adapting the three- 
dimensional expressionsu to two dimensions we have 
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M^^ is the inverse mass tensor with the elements 
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V the slab volume, w„(fc) — ■|Ve„(fc) the electron group 
velocity, e„(fc) the energy of an electron with momentum 
k in the subband ri, and / the Fermi function. The sums 
in Eqs. (Q) and (^) run over the wavevectors k in the 
two-dimensional hexagonal Brillouin zone and over the 
electron subband index n. (Each atomic electron level 
splits into d two-dimensional electron subbands where d 
is the number of monolayers in the slab.) The spin degen- 
eracy is included in the prefactors. As will be seen below, 
the relaxation time Tn(k) of thin slabs at low tempera- 
tures is dominated by elastic scattering on the thickness 
fluctuations of the Pb film (roughness of the interface to 
the substrate). Therefore, we will take the limit of short- 
wavelength surface height fluctuations (with their corre- 
lation length ^ ^ 1/fcF, kp being the Fermi wavevecij 
tor). In this case the relaxation time is fc— independent .□ 
Since no detailed information on the scattering mecha- 
nisms and on the interface structure is available, we as- 
sume that T is also independent of the subband index n. 
Then, Rh is independent of r and is determined solely 
by the electron bandstructure at the Fermi energy ep. 

The electronic bandstructure calculations were per- 
formed on relaxed 1 to 8 monolayers thick (111) oriented 
Pb slabs with the in-plane lattice constant a — 3.50 A, 
separated by « 10 A of vacuum. The electron band 
energies and the total energy have been calculated us- 
ing the full-potential linearized augmented plane-wave 
method within_Lhe WIEN2k codalj in the local-density 
approximation.^ A Pb muffin-tin radius of 2.6 a.u. and a 
tetrahcdral mesh of 623 /c— points in the irreducible part 
of the full Brillouin zone were used in the self-consistent 
electronic structure calculations. The kinetic-energy cut- 
off was set to E^^^ — 7.3 Ry and the plane-wave expan- 
sion cutoff to — I^y- Later, in the magnetocon- 
duction calculations, we used about 4000 k points in the 
irreducible part of the two-dimensional Brillouin zone. 
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FIG. 1: Fermi lines of relaxed 3 and 6 ML thick free-standing 
Pb slabs. The solid lines encircle electron and dashed lines 
hole pockets. 



Several electron subbands cross the Fermi energy and 
there are several Fermi lines in the two-dimensional Bril- 
louin zone. As typical examples. Fig. ^ shows the Fermi 
lines of 3 and 6 ML thick slabs. The number of Fermi 
lines is an increasing function of the slab thickness. Only 
electrons close to the Fermi lines take part in the trans- 
port and cto/t and <jh/t'^ are obtained after numerical 
integration of expressions (||) and (^) along the Fermi 
lines. Whereas ctq is always proportional to the density 
of states at the Fermi lines , au of a subband is negative 
if the Fermi line makes a closed loop around occupied 
states (electrons) and positive if the loop encircles empty 
states (holes) in a subband. It is the balance between the 
number of electron and hole states at together with 
their effective masses which determines the sign of Rh- 
As a consequence, Rh is sensitive to the details of the 
bandstructure. 

Fig. |2| shows the experimental conductivities cr™^ of 
annealed Pb films on Si(lll)— (7 x 7) together with the 
data of epitaxial Pb on Si(lll) covered with Ag.Q Three 
properties of these data indicate that interface and/or 
surface scattering is the dominant mechanism of (elastic) 
electron scattering in this system. Firstly, only a weak 
temperature dependence of a"^^ is found below 100 K. 
Secondly, there is only a weak dependence of conduc- 
tivity on crystallinity and order in the films. Anneal- 
ing of amorphous films grown at 15 K resulted in an 
increase of conductivity, but this increase never exceeded 
a factor of 2.5. Thirdly, the roughly linear thickness de- 
pendence of cr™^ is again consistent with the assump- 
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FIG. 2: Experimental conductivity ao and the calculated av- 
erage relaxation time r. The symbols are for annealed Pb 
films on Si(lll) — (7 x 7), measured at three different tem- 
peratures, the dashed line is for eMtaxial Pb on Si(lll)- 
Ag(V3 X ^/3) at T = 20 K from RefB and the solid line is a 
linear fit to the low-temperature data for d < 6ML. The main 
process limiting the conductivity and the relaxation time in 
Pb ultrathin films is scattering on the interface fluctuations, 
the disorder within the film plays only a minor role. 



tion just made that the dominant scattering mechanism 
in these uhr|ai-tliin films is diffuse scattering at interface 
fluctuations.BEJ For thicker samples the scattering at im- 
perfections within the film becomes noticeable, and the 
conductivity starts to deviate from the linear behavior 
and eventually saturates in the bulk value, where the 
scattering on lattice defects and at higher temperatures 
also on phonons limit the conductivity. 

We measured aQ^^ for arbitrary film thickness d. At 
non-integer thicknesses the film surface is covered with 
islands and their borders (i.e., steps) act as scattering 
centers for the electrons. Therefore one would expect 
a dip in aQ^^ for non-integer thicknesses. However, a 
small effect of this kind has only been observed at c? > 12 
ML. This suggests that the typical island sizes exceed 
the charge-carrier mean free path for very thin layers 
and electron scattering at the vacuum interface is not 
important for films with d < 10. 

Clearly seen is an offset in the conductivity towards 
d ~ 1 ML. The onset of conductivity due to percola- 
tion of the evaporated Pb at a surface temperature of 15 
K was found to be at 0.8 ML. The very low conductiv- 
ity of this first monolayer indicates that the main trans- 
port mechanism in this layer is hopping between localized 
states. Indeed, the conductivitw shows a temperature- 
activated behavior for d — I MLB The first Pb monolayer 
on Si(lll) — (7 X 7) has the periodicity of the substratecl 
and may incorporate the adatoms of the (7x 7) substrate. 
This might be the reason for the very low conductivity 
and for the localization of charge carriers. 

The close agreement of the conductivity data obtained 
in our measurements with those of Ref. |^ is not at all self- 



evident, since our measurements have been carried out on 
Pb films that are expected to be much more disordered 
than those grown epitaxially on Si(lll)-Ag(-\/3 x V3)- 
Since the thickness dependence is almost the same for 
Pb films prepared on both substrates, we must conclude 
that interface scattering plays the same dominant role on 
both substrates. This means on the other hand that the 
degree of order in these ultrathin films is of secondary 
importance for electronic transport. We also note that 
the Ag interface layer obviously does not contribute to 
the measured conductance. 

In contrast to the experimentally measured (T™^, the 
calculations of ctq /''' and an / were performed only for 
integer d. Surprisingly, cfq/t turned out to be almost in- 
dependent of the slab thickness in the investigated range 
between 2 and 8 ML, 
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With increasing number of monolayers the number of 
subbands increases and thus the conductance but not the 
ratio (Tq/t. This means that the contributions of newly 
emerging subbands to CTo /''' as a function of layer thick- 
ness are very similar. Eq. (^) implies that r is propor- 
tional to Co- Using the experimental values for ao, we 
determined t, also shown in Fig. ^j. 

It is interesting to look at the charge-carrier mean free 
path A = vpT. From the bandstructure we get the Fermi 
velocity vp ^ 1.2 x 10^ m/s and, using the values for 
T from Fig. |[ we see that A is of the order 1 to 2 film 
thicknesses. It is therefore tempting to say that the elec- 
trons scatter on the film surfaces. However, the charge 
carriers scatter only on imperfections (steps, misplaced 
atoms etc.) on the surfaces and not on a perfectly flat 
surface. For d = 1 ML the mean free path is shorter than 
the lattice spacing and the wave- vector description of the 
electron states breaks down. 

The relaxation time causes broadening of the electron 
energy levels in the vicinity of the Fermi energy. There- 
fore, Rh was calculated as {ao/r) / {an /t"^) where () is 
an average over the Gaussian distribution with the half- 
width S ~ h/r. 

The experimental results (open symbols) of the Hall 
coefficient are compared in Fig. ^ with those calculated 
(solid circles). Whereas the calculations have been done 
for perfectly periodic slabs without substrate, the mea- 
surements refer to epitaxial films with different degrees of 
disorder. Clearly seen are the anomalous oscillations of 
Rh as a function of film thickness d for both sets of exper- 
imental data. Our experimental results (squares) agree 
very well within the experimental error with the results 
of Ref. ^ (open circles in Fig. ||). After having found 
close agreement for the data of electrical conductance, 
the agreement of the results for Rh are not surprising, 
since the disorder within the film plays a negligible role. 
On the other hand, this result justifies our assumption of 
scattering times independent of k and n. 

Rh, calculated in this paper without any adjustable 
parameter, reproduces the oscillating behavior seen in 
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FIG. 3; Quantum-size induced oscillations of the Hall coef- 
ficient with the thickness d. Solid circles: calculations for a 
free-standing Pb slab; open symbols: experimental values for 
Pb films on two different Si(lll) substrates (open circles are 
from Ref. ^ Unless explicitly indicated, the experimental er- 
ror bars are of the size of the symbols. The estimated error 
bars of calculated Rh are ±1 x 10~^^ m'^/C and come mainly 
from uncertainty in the broadening of the energy levels due 
to scattering. The solid line is a guide for the eye through the 
experimental points. 

the experiment semi-quantitatively and demonstrates 
that the main reason for this behavior are the thickness 
dependent changes in the band structure. The devia- 
tions at d — 2 and d — 5 may be a consequence of the 
approximation of equal relaxation times for all subbands. 
Especially in the case d — 2 also the 7x7 periodicity of 
the substrate, which was neglected in the slab calcula- 
tion, might play a role. Our calculations also show that 
fJo /t is a slowly varying function of energy. Therefore it 
is not very sensitive to the exact location of the Fermi 
energy within the band structure, and the conductance 
increases smoothly with film thickness, and with only 



small variations during deposition. In contrast, the cal- 
culated cth/t^ is a strongly energy-dependent function 
and most of the energy dependence of Rh comes from 
(Jh/t'^, therefore Rh is sensitive to averaging. 

To conclude, we have shown that only a quantita- 
tive bandstructure analysis, done here with DFT, al- 
lows an understanding of the anomalous magnetocon- 
ductive properties of ultrathin Pb films. Oscillations in 
Rh as a function of film thickness are the consequence 
of competing - and to a large extent compensating - 
contributions of electrons from different subbands and 
cannot be explained with any model which does not 
take into account the details of the electronic bandstruc- 
ture. The relaxation time and thus the interface fluctu- 
ations play a subordinate role, they affect Rh through 
the width of the distribution function and through the 
subband-dependence of the relaxation time, which was 
neglected in this paper. Although our calculations with 
free-standing films already show all qualitative features 
observed in experiment, it is obvious that this method 
could be further refined by explicitly taking into account 
the underlying substrate. This might help to understand 
the specific role of the metal-semiconductor interface con- 
sidered here. Generally, this interplay between calcula- 
tions and experiments is of primary importance since in- 
ner interfaces cannot be varied easily in experiments. In 
this context, also the open problem of the very low con- 
ductivity and the high Rh of the first monolayer may be 
solvable. This involves most likely hopping conduction 
between localized states in the Pb film on reconstructed 
Si(lll) substrate. 
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